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NOMENCLATURE 


Surface area of cylinder exposed to free stream air, те. 
у іпбегсері. 
Multiplication constant, dimensionless. 
Diameter of cylindrical element, m. 


Roughness Reynolds number based on e, a characteristic 
height of the roughness, dimensionless. 


Heat transfer similarity function, dimensionless. 
Local heat-transfer coefficient, W/m?-?C. 
Distance between channel splitter plates m. 
Current drawn, A. 

Thermal conductivity of free stream air, W/m-?C. 
Height of cylindrical element, m. 

Slope of a line. 

Row position number in array. 

Nusselt number, hH/k. 

Prandtl] number, Cpu/K, dimensionless. 
Power input from electrical power supply, М. 
Power lost from cylindrical elements, W. 
HU o / у 


Reynolds number, , dimensionless. 


Surface temperature of the cylinder heated by the 
resistor, °С. 


Surface temperature of the ith cylinder heated by the 
thermal wake, °С. 


D Freestream air temperature, “C. 
U Freestream air velocity, m/s. 

V Power supply voltage, V. 

Greek Symbols 


9 Dimensionless temperature of a cylinder in the same . 
column as the cylinder heated by a resistor. 


V Kinematic viscosity, m@/s. 
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I. INTRODUCTION 


The removal of heat from electronic equipment to assure 
reliable operation has been a problem plaquinq the electronic 
industry for many years. As a result, many researchers have - 
attempted to develop specific correlations to enable designers to 
attempt to predict possible solutions to this problem. 

This research concerns air flowinq across cylindrical 
elements mounted in an array on one wall of a vertical adiabatic 
channel. This configuration is much the same as an array of 
transistors in the T05 configuration arranged on stacks of 


parallel circuit boards in various types of electronic equipment. 


A. THEORY OF VERTICAL CHANNELS FORMED BY CIRCUIT BOARDS 

Many electronic systems consist partially of arrays of 
printed circuit boards supporting a fairly dense array of 
electronic elements on one side as shown in Figure l.1. These 
circuit boards can be arranged so that the side with the elements 
(front side) faces the relatively smooth backside of another. 
When a series of circuit boards are arranged in such a fashion 
(see Figure 1.2), long vertical channels in which cooling air can 
be supplied by forced convection can result. During operation, 
the elements such as flatpacks, dips resistors and transistors 
generate heat, the majority of which is dissipated to the 
surroundings with little transmitted back through the circuit 


board. As a result, it appears as if all of the heat is 


121 





Figure 1.1 Modern Circuit Board 
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Modular Arrangement of Circuit Board 
13 


с лге 1.2 


transmitted by one side of the channel or one side of each 
Circuit board to the air stream or air layer between the boards. 

Complexity is added to the system when elements serve 
different purposes and, as a result, possess different sizes, 
shapes and heat dissipating characteristics. 

The process of heat removal from electronic systems creates a 
very comlex design problem governed by such factors as circuit 
board/element arrangement, spacing, heating characteristics and 
operating environment. Moreover, the ultimate goal of any 
thermal analysis and the prediction of operating temperatures is 
the assessment of reliability. Reliability is the probability 
that an electronic component/system “111 operate properly without 
failure for its intended lifetime under its environmental 
operating conditions. [Ref. 1] 

For many years, there have been numerous attempts to derive 
techniques to aid designers in this field. Unfortunately, the 
methods developed are often too complex, inadequate or restricted 
to specific condi tions. 

Early work with natural convection between parallel arrays of 
vertical circuit boards was performed by Elenbaas [Ref. 2], 
Bodoia and Osterle [Ref. 3] and Engel and Muller [Ref. 4]. In 
1972, Aung et al. [Ref. 5] studied the effects of circuit board 
supports and elements in the field of flow for wide channel 
Spacings. Aung went on to investigate uniform heat flux on two 
channel walls in which he subsequently studied the effects of 


baffles and circuit board staggering. [Ref. 6] There is 
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considerable literature describing the efforts that have been 
made to attempt to model circuit channel cooling by natural 
convection and forced convection. 

For example, Kelleher utilized a thermal network approach to 
derive a more representative model of forced convection for one 
specific geometry with equal heat releases [Ref. 7). Other 
experimental or semi-experimental approaches were performed by 
Tamura [Ref. 8], Finch [Ref. 9], Laermer [Ref. 10], Marto [Ref. 
11] and Parkes and Preston [Ref. 12]. Most of this work is 


mestricted to specific conditions. 


B. THEORY OF SURFACE ROUGHNESS GEOMETRIES 

Surface roughnesses have been classified by several different 
techniques of which the most commonly used are the k and d type 
roughness classifications consisting of 2 and 3 dimension classes 
which are further divided into artificial and commercial 
roughness subclasses. Although many researchers allude to these 
classifications, only a limited number actually use this 
terminology in agreement. It is further noted that in the area 
of heat-transfer, the k- and d- type roughnesses seem irrelevant 
because the same heat transfer roughness parameter (g*) 


correlation is proposed for all types of roughnesses. 


C. THEORY OF HEAT TRANSFER FROM ROUGH SURFACES 


Earlier researchers in this field seem to agree that the 


heat-transfer roughness parameter (g*), commonly referred to as 
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the inverse of the sublayer Stanton number, takes the functional 
form: 

a Н р q 

g = l/Stsublayer - C(e ) (Pr) 
where e* is the roughness Reynolds number based on an arbitrary 
characteristic height of the roughness (e). The different values 
of p and q fall into two categories: 

1. The assumption of similarity near the wall which results 
in value of p = .2 and q = .44 or for the eddy-diffusivity 
models p = .25 and q = .66 [Refs. 13 and 14]. 

2. The assumption of a breakdown in close to the wall 
similarity results in values of p = .45 and q = .57 or 
р = .5 and q = .66 for the eddy diffusivity models 
[Refs. 15 and 14]. 

With regard to this study, the foregoing is only of academic 


interest because the geometry here is defined as an in-line array 


of cylindrical elements. 
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II. EQUIPMENT DESIGN AND INSTRUMENTATION 


The equipment utilized in this research was either 
manufactured by the Engineering Machine Shop or supplied by the 
Naval Engineering Department of the Naval Postgraduate School, 


Monterey, California. 


A. WIND TUNNEL 

The wind tunnel (Figure 2.1) is constructed of 1/8" 
plexiglass with the exception of two 1/4" access plates and 
various reinforcement members mounted on the exterior so as to 
ensure smooth air flow on the tunnel interior. 

The test section of the tunnel has a 12" by 4" cross section 
and a height of 18" with access covers on both faces to 
facilitate instrumentation and support for the test surfaces. 

The tunnel intake is belled with an 18" 201107 to produce an 
orifice that has a 40" by 12" cross section. An adjustable leq 
is attached at each corner to serve as a support for the tunnel, 
to level the tunnel and to provide an 8" clearance off the floor. 

| The upper end of the tunnel is reduced to a section of 2" PVC 
ducting which is connected to a 550 cfm turbo fan. Air flow 
throuah the tunnel is controlled by a series of valves connected 
in line with the fan. 

A 1/8" plexiglass splitter plate was mounted 1/2" off of the 
access plate facing the test surface in order to reduce the 


effects of the boundary layer formed in the tunnel. This 
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distance was pre-determined through the use of computer software 
[Ref. 16] and the results are displayed in Figure 2.2. A wooden 
frame box (4'x4'x 8') was constructed around the wind tunnel to 
elliminate the effects of incident radiation on the cylinders 
from the laboratory environment and the box also served to reduce 
any transient air currents created by movement in the room (i.e., 


doors opening and closing). 


BE TEST SURFACE 

The design of the test surface used in this research was 
derived from the previous work of Arvizu and Moffat [Ref. 17] and 
Ortega [Ref. 18]. It was constructed of a 12" x 12.5" x 1/8" 
sheet of plexiglass (for support) with a 12" x 12.5" x 3/8" sheet 
of balsa wood (thermal boundary) cemented to it. The leading 
edge of this (balsa wood-plexiglass) splitter plate was cut at a 
35° angle to form a sharp leading edge to reduce the effects of 
leading edge separation. The surface (balsa wood) of the 
splitter plate was then sanded and painted three times with white 
enamel paint and sanded a fourth time to ensure a smooth test 
surface. Small (1/16" diameter) holes were drilled through the 
plate to permit the installation of resistor leads and 
thermocouples into the cylindrical elements mounted on the balsa 
wood. Three mounting screws were attached to the plexiglass side 
of the splitter plate and to the access plate to facilitate the 
positioning of the test surface. A tail fin mounted to the tail 


end of the splitter plate was used to regulate the air flow on 
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both side of the splitter plate. A fourth adjusting screw was 
attached to this fin and through the access plate 

(Figure 2.3). 

Р CYLINDRICAL ELEMENTS 

The cylindrical elements used in this research were 
manufactured from 3/4" round aluminum stock and cut to a height 
of 1.0 cm. A 7/16" diameter hole was cut to a depth of 1/4" in 
the center of the element and two additional holes (1/8" 
diameter, 1/4 deep) were cut opposite each other to permit the 
installation of the 56 ohm, 1/2 watt resistors. Another hole 
(1/32" diameter, 3/8" deep) was drilled to position the 
thermocouple (Figure 2.4). | 

The resistors were epoxied into forty two of the cylinders 
and six more were cut in half so that an array of /x6 cylinders 
could be mounted evenly on the splitter plate and the 12 half 
cylinders could be arranged along the edges of the plates (Figure 
21). 

The resistor (56 ohm, 1/2 watt), prior to insertion into the 
cylinders, had the leads cut and smaller leads (2841/1, #30 solid 
spc.) attached so as to reduce heat losses by conduction through 
them. 

The cylinders were mounted on the test surface in an in-line 
array with a separation of one cylinder diameter between 


cylinders. This was done both vertically and horizontally. 
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Top Row: Assorted Views of Cylinders 
Bottom Row: TO5 Transistor and a nickel 


Figure 2.4 Photograph of Cylinders. 


T) 
j^. 
5 





Figure 2.5 Photograph of Cylinders Arranged in an Array 
оп the теза тасе 


D. INSTRUMENTATION 
1. Temperature Measurements 

Type T, copper-constantan thermocouples were utilized for 
the determination of the temperature of the heated elements and 
the inlet and outlet temperatures of the wind tunnel near the 
test surface. The thermocouples were connected to seven scanner 
boards in an autodata nine and calibrated to an accuracy of 
+0.2°F (Appendix A.) The thermocouples were then routed through 
an access plate in the tunnel through the splitter plate and into 
each of the elements. 

2. Velocity Measurements 

The velocity measurements were obtained by the use of a 
micro manometer and a pitot-static tube. Changes in the dynamic 
and static pressure were measured by the micro manometer and 
using Barnoulli's equation, velocity was determined. The pitot- 
static tube was placed between the splitter plate and the test 
surface splitter plate and used to determine the approach 
velocity between them. 

A TSI (Model 1050/1050AA) hot wire anemometer was used to 
determine the velocity profile prior to entry into the test area. 
The hot wire anemometer was calibrated with the micro manometer 
and used only to ensure that the approach velocity was uniform. 

Accuracy of the mirco manometer was determined to be 
*.0005 inches of water and resulted in an error of *0.2 ft/sec or 


+7Ст/$. Because the micro manometer was utilized for the 


Zo 


calibration of the hot wire anemometer, more accuracy might 
result if the micro manometer was used to determine the velocity 
between the splitter plates. 

3. Power Measurements and Control 

A LAMBDA (60 volt, 10 amp) Regulated Power Supply (Mode! 
LX 345AFM) was used to power to the 56 ohm, 1/2 watt resistors 
mounted in the forty two cylindrical elements. 

A switch board was set up to allow power to be delivered 
to individual elements in either of two columns, the third or the 
fourth columns. It was also set up in such a manner that each or 
all rows could be energized. 

A Hewlett Packard digital multimeter (3466A) was used to 
measure the voltage drop across the resistors and a Fluke 80604 
multimeter was used to measure the current through the circuit. 
The power delivered is the product of current and voltage. 

4. Data Acquisition 

The Autodata Nine mentioned above is a self contained 
recording device with up to 100 channels available for both 
temperature and voltage measurements. During this research, 71 
channels were utilized, one for voltage measurement and the rest 


for temperature measurements in degrees Fahrenheit. (Figure 2.6) 


E. EXPERIMENTAL PROCEDURE 
For this research, data was recorded at three different 
channel height to cylinder height ratios (1.0, 2.3 and 4.6) for 


two different approach velocities (1.5 m/s and 3.0 m/s). The 
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Figure 2.6 Photograph of Autodata Nine. 
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data recorded consisted of the voltage supplied to the circuit, 
the current through the circuit and the temperatures of the test 
elements as well as the inlet and exhaust temperatures of the 
airstream. For each channel height and approach velocity, data 
was gathered after a steady state condition was reached, for . 
each element of a column energized one at a time and for the 


first four rows and for all six rows energized all at once. 
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III. DATA REDUCTION 


A. DATA TAKEN 
Net bree Stream Temperature, T (°F) 
Cylinder Surface temperature, Toyo (°F) 
Atmospheric Pressure, Patm [inches Hg.) 
Pitot-Static Pressure Difference, AP (inches Н20) 
Power Supply Voltage, V (volts) 
Current Drawn, I (amps) 


uS CF) 


B. SOUGHT AFTER PARAMETERS 
l. Dimensionless Temperature Ratio 
The Dimensionless Temperature Ratio, Ө, is obtained as a 
function of cylinder location, channel configuration and 


freestream velocity. Here 6 is defined as 


6 A Tei сащ 


аж. Доц Ta 
where Ts,1 is the surface temperature obtained by the first 
passive cylinder directly after the cylinder heated by a 
resistor. Те,і is the surface temperature obtained by individual 
passive cylinders after and in the same column as the cylinder 
heated by the resistor where i is the row after the heated 
cylinder. Note that 
MS 


5.17 Те 


- 1.0 by definition 
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2. Heat-transfer Coefficient 
The Heat-transfer Coefficient, н(м/т с), is found and 
is a function of the cylinder location, channel configuration and 


freestream air velocity. 
The heat transfer coefficient was determined as follows: 


(W) - с (W) 
h = input losses 
А (1560 - T. 

Ts,o and T, were recorded by the autodata nine and 
corrected by the calibration data in Appendix A. 

A, is the surface area of the cylinder is calculated for 
channel/cylinder height ratios of 1.0, 2.3 and 4.6. If the 
channel/cylinder height ratio is equal to 1.0, A = mdL and if the 
channel/cylinder height ratio is greater than 1.0, А * m[dL + 
42/4]. 

The value of Q159sseg in watts is derived from a parabolic 
temperature dependent approximation utilized by Ortega [Ref. 18] 
in his researach with aluminum cubes. This approximation is 


assumed to be the same for cylinders and is used here. 


Oosses (М) з 20962927 X103 000950599 00212778 
- 0.023695 x 10-5 (T,,g-T, )? 


Ortega also determined the radiation losses to be an order of 
magnitude smaller than those given above, and radiation losses 


were therefore disregarded. 
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Qinput(W) is the product of power supply voltage and 


current drawn. 
Qinput = V(volts) x I(amps) 


3. Reynolds Number 
The Reynolds number is developed as a function of channel 
configuration and free stream air velocity. In this research, 
the Reynolds number was calculated using the freestream air 
velocity, channel height and viscosity as follows: 


ет = Channel height x velocity 


viscos1ty 


The value of the viscosity is obtained by a linear 
interpolation of data from tables of viscosity and air 
temperature. 

4. Nusselt Number 

The Nusselt number is prescribed as a function of channel 
configuration. Here, the Nusselt number was calculated using 
channel height, heat-transfer coefficient and the thermal 


conductivity of the free stream air as follows: 


uy - channel height x heat-transfer coefficient 


thermal conductivity 


The value of the thermal conductivity was obtained by a 
linear interpolation of data from tables of thermal conductivity 


and air temperature. 
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C. TYPICAL БЕТ 28 227. 
Table I presents a typical set of data obtained. These data 


will be used here to illustrate the data reduction procedure. 


TABLE I 
DATA SET 


Patm = 29.736 inches Hg. 

AP = 0.0051 inches H>0 
Channel/Cylinder Height Ratio = 4.6 
Cylinder Heated - row 1, column 3 


Cylinder Surface Temperature, T.,g = 206.1°F 


Voltage Supplied 14.380 volts. 


Current Drawn = 0.225 amps. 
Tes ORs 
Te AE COND 
Теза = ишао r 
Тед ~ S F 
Less > и 
T = OCF 


со 


D. CALCULATIONS 
1. Dimensionless Temperature Ratio, 
The temperatures for the cylinder surface were recorded 
by the autodata nine and corrected by use of the thermocouple 


calibration data presented in Appendix A. 


Тс, (corrected) = M X Iç,; (recorded) © 
N RECO CET b Ts,i(Corrected) 


Ts; = 1.00267909073) = 0.908593 = о 
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ОЕЕО ЕО - 0.19287 = 73.8 Е 
EN UE (7287) - 0405625 s 72.3 F 
Ts»4 = 1,0018 (68.6) - 0.5980 = 68.1 °F 
Des Tas 
A 
Toes] ын Tes 
01 = 79.6 - 68.1 = 11.5 = 1.0 (by definition) 
P eos. 5: 656 
79.6 € 68.1 11.5 
95 з 929 2277 = 0.3652 
о = • 11.5 
2 - 68.1 = 3.1 = 0.2696 
= 70.3 - 68.1 = 2.2 = 0.1913 
79.6 = 68.1 11.5 


2. Heat-Tansfer Coefficient 


поз Qinput - Qlosses 

A с» () = со 
ШО (010025 (206.1 Р) - 0.4745 
Т = 1.0018 (68.6 F) - 0.5980 


со 


Шт. — 206.1 = 658.1 = 138.0 


Channel/cylinder height ratio = 4.6 


F 


206.1 Е (corrected) 
BE IF (corrected) 


= 76.7 С 


therefore A = п[(0.019т) (0.01т) + (0.019т) 2/41 
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А = 0.00088т2 
Отоѕѕеѕ = 046927 х 10 s ОШ ИШ и 
- 0.23695 х 10-601..0-Т.)2 


where Tso -Iœ is iN "C and Qlosses 15 in watts. 


Те. о“ Гоу = 7168896 

n - 0.373 watts 

ш = VxI = 14.380 volts x 0.225 amps = 3.27 watts 
h = 3.27 - 583. - 42.92 W/m“ °C 


3. Reynolds Number 


Rey = HU./ Y 


H = 4.6 (.01т) = 0.046m 


Ше 1.5m/s 


38,02 х 10-6 = 109.9 х Оа гт? 


е е o 
where Ша is in С 


T = 20.0°C 


38.02 x 10-9m?/s 
Rey = .046m x 1.5m/s = 1814.7 
38.02 x 10-6m2/s 


4. Nusselt Number 


Nuy = Hh/k 
Н = .046m п = 42.92 W/m?'C 
К = 0.0407 - 76.8 х 10-6 (27.0-Т ) 
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. . o 
uere T. is in С 


k = 40.162 x 10-3w/m°C 


Миң = .046m x 42.92 W/m2 C = 49.16 
40.162 x 10-3 W/m°C 
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IV. RESULTS 


The experimental measurements that were obtained lead to two 
different subject categories: thermal wake effects and 
heat-transfer coefficient. Table II presents the range of the 
parameters covered and a listing of the data is tabulated in 


Appendix C. 


TABLE II 
PARAMETERS TESTED 


Thermal Wake 
Channel Height/Cylinder Height Ratio l. 
Downstream Row Number О е 
Velocity (m/s) 15558 
Heat-Transfer Coefficient 
Channel Height/Cylinder Height Ratio 1l 
Heated Cylinder Row Number 1.2 
Number of Rows Heated 4 
Velocity (m/s) 12:52 
А. THERMAL WAKE 

The thermal wake created by a single heated cylinder in an 
array is influenced by the channel height to cylinder height 
ratio, the mass flow rate and the position of the cylinder in 
the array. 

For a channel height/cylinder height ratio of 1.9, there, is 


no free stream, and the flow through the array undergoes periodic 


accelerations and decelerations. It therefore seems unlikely 
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mat da channel Wwerght/cylinder height ratio of 4.6 would possess 
the same convective means for determination of the downstream 
wake. However, a plot of & as a function of downstream row 
number (Figure 4.1 and 4.2) demonstrates a 1/N dependence for 
both velocities at both channel heights. 

Figure 4.3 shows an even closer resemblance to the 1/N curve 
for the thermal wake created by a heated cylinder located in a 
row internal to the matrix of cylinders. 

The following is provided for the understanding of the 
graphical curves found in the legends of various graphs in 


Chapter IV: 


Hx .XRxVx.x 

Hx. x Channel Height/Cylinder Height Ratio x.x 

Rx Row that the heated cylinder is in x 

X X Approach Velocity x.x 
For example, H.4.6R1V3.0 denotes data gathered at a channel 
height/cylinder height ratio of 4.6, where the cylinder heated 


is in the first row and the approach velocity is 3.0 m/s. 


B. HEAT-TRANSFER COEFFICIENT 

In Figure 4.4 the heat-transfer coefficient, h, for a single 
heated cylinder in the first row is plotted against channel 
height/cylinder height ratios for approach velocities of 1.5 m/s 
and 3.0 m/s. It can be seen that as the channel height 
increases, the value of h drops off fairly rapidly to an 


asymototic limit at a channel height/cylinder height ratio of 
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4.6. Air flows over the cylinders instead of through them (the 
path of least resistance). 

Figure 4.5 show a plot of h against the channel height to 
cylinder height ratio for an internal cylinder where h drops off 
a little more rapidly. 

Figure 4.6 plots h against array position for a channel 
height/cylinder height ratio of 1.0. Here, h increases with 
position upstream. This correlates to the increased turbulence 
with position downstream. 

Figure 4.7 is a plot of h against array position for a 
channel height to cylinder height ratio of 4.6. For this 
configuration h appears to increase slightly and then drop off 
Slightly so that the value of h at the sixth row is less then 
that of the first row. 

Figure 4.8 plots Nusselt number against Reynolds number for 
all three channel height to cylinder height ratios. This 
composition leads to a correlation of the form Nu = C Re0.55 
where C is different for different array densities as noted by 


Arvizu and Mof fat [Ref. 17]. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


In this research, for cylindrical elements mounted in a 
vertical channel, the characteristic curves produced closely 
resemble those for a similar configuration of cubical eI ements 
which were investigated by of Arvizu and Mof fat [Ref. 17]. 

Although there is room for more data to assist in making firm 
conclusions about the flow characteristics of air flow around 
cylinders, the 

The following has been observed: 


1. The 1/N characteristic curve produced by the thermal wake 
of a heated cylinder on an adiabatic wall in a channel, 
may allow for fairly accurate prediction of temperature 
increase in downstream elements in an array for various air 
flow rates. 


2. Channel spacings and airflow rates have very little affect 
on the shape of the characteristic curve produced by the 
thermal wake of a heated cylinder in an array. 


3. The position of a heated cylinder in an array, on an 
adiabatic wall in a channel, has a definite affect on its 
heat transfer characteristics. Elements requiring more 
cooling should be placed closer to the freestream air inlet. 
This would reduce effects created by the thermal wake of 
other cylinders. 


4. Heat transfer in an array of cylindrical elements is 
mainly influenced by the air flow between the cylinders 
and slightly less by the turbulence in the air flow. 

5. As the channel height to cylinder height increases, the 
amount of heat transfer that occurs decreases to an 
asymototic limit near a channel/cylinder height ratio of 
4.6. 

For continued work in this field, it is recommended that, a 


longer test surface with more rows of elements be utilized along 
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with higher velocities and various geometric shapes. All of this 
can be accomplished with only a few modifications to the 


equipment at hand. 
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APPENDIX A. 
THERMOCOUPLE CALIBRATION 

The thermocouples were all connected to a wooden framed 
Switchboard which was connected by 10' wire leads to seven 
scanner boards mounted in the autodata nine data acquisition 
device. This combination was then calibrated as a complete unit 
in the calibration laboratory against five known temperatures: 

66.02°F, 106.14°F, 127.615°F, 151.20°F and 200.42°F 


Thermocouple readings follow: 


35 66.3 1065 127.9 151.2 200.4 
CD Tw EVO 75. 704 36 66.3 106.5 | 128.0 | 151.4 | 200.5 
: а До S cor сор 37 66.3 106.5 127.9 151.3 200.5 
По 1220 151. alb зв 66.4 106.5 128.2 181.3 200.6 
à Mc 2123. Sa e 39 566.3 106.6 127.9 | 151.2 200.4 
і цог 277 22212: 222 40 66.6 106.7 128.2 151.6 200.7 
м DOM 7 нг 41 66.5 106.6 128.1 151.5 200.7 
: СО lena 15:41: "V 42 66.6 106.8 128.4 151.6 200.9 
: Dr db ls bio 208 43 66.6 106.7 128.2 1515 200.8 
2 22 eg з 127-700 510 75234 46 66.7 106.8 128.4 151.6 200.8 
10 66.2 106.6 127.8 151.2 200.4 2 МИ 222272: 18021 
11 66.8 106.9 128.5 151.8 201.0 И A кос в 
12 66.8 107.0 128.5 151.8 200.9 "EM uc bez сеу 2:098 
13 66.7 106.8 128.3 151.6 200.8 " E Noir WEE. лан соб 
14 66.7 107.0 128.6 151.8 201.0 E Огоо Мао зїг собі 
15 66.6 106.8 128.3 151.6 200.7 5 E eu docs М e 
16 66.7 106.8 128.3 151.7 200.8 - т о A 
17 66.5 106.7 128.0 151.4 200.7 a IN аа с а ШЕК, 
18 66.7 106.8 128.2 151.6 200.8 2 "NEUEM M. 
19 64.5 106.6 128.0 151.4 200.7 2 Mo м №. 
20 64.2 106.8 122.6 182.6 199.3 Ё о mM UM OM - 
21 66.2 106.6 128.0 181.3 200.4 E Mu LU РВЕ 
я з e: 120 = 20005 57 66.6 106.5 128.0 151.6 200.6 
НА гал oo Ro. 200.6 S8 66.6 106.6 128.2 151.4 200.6 
МИН О o БОА 200.5 59 66.5 106.6 128.1 151.6 200.6 
ИГ пр Дейв gene. 200.4 60 — 66.2 106.8 128.0 151.5 200.4 
SS 19Щес с 2906 61 66.2 106.7 128.0 151.4 200.5 
ET A e 127-1 151.0 200.2 62 66.3 107.1 128.2 151.4 200.6 
ПИ се A ssp 200 63 66.3 106.8 128.0 151.5 200.6 
ВИ М с: 1272 151.1. 200.3 4 66.5 106.9 128.20 151.4 200.6 
= IEEE S О 20041 55 66.3 106.8 128.1 151.4 200.6 
Р} Банн 1214 ous 2664 66 66.3 106.8 128.1 151.4 200.6 
SEE oec uem 290:6 67 66.4 106.9 128.5 151.5 200.7 
27 -—— OL NE ME. 88 66.3 106.8 128.1 151.5 200.7 
34 66.3 106.5 128.0 151.4 200.6 
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A least squares fit program was developed to produce 


thermocouple corrections with the following result: 


Tactual 
THERMOCOUPLE 1 мъ 
THERMOCOUPLE 2 Ма 
THERMOCOUPLE 5 мъ 
THERMOCOUPLE 4 Ma 
THERMOCOUPLE 5 Ma 
THERMOCOUPLE 6 ма 
THERMOCOUPLE 7 ма 
THERMOCOUPLE 8 Ма 
THERMOCOUPLE 9 ма 
THERMOCOUPLE 10 Ms 
THERMOCOUPLE 11 ме 
THERMOCOUPLE 12 на 
THERMOCOUPLE 13 Ma 
THERMOCOUPLE 14 Ma 
THERMOCOUPLE 15 Ma 
THERMOCOUPLE 16 Ms 
THERMOCOUPLE 17 Ms 
THERMOCOUPLE 1B Mu 
THERMOCOUPLE 19 Ha 
THERMOCOUPLE 20 Ms 
THERMOCQUPLE 21 Ms 
THERMOCOUPLE 22 Ms 
THERMOCQUPLE 23 Ma 
THERMOCOUPLE 24 Ma 
THERMOCQUPLE 25 мэ 
THERMOCOUPLE 26 Mm 
THERMOCOUPLE 27 Мэ 
THERMOCOUPLE 2B Мэ 
THERMOCOUPLE 29 ча 
THERMOCOUPLE 30 Ms 
THERMOCOUPLE 31 Ма 
THERMOCOUPLE 32 мъ 
THERMQCQUPLE 55 Ма 

34 М» 


THERMOCQUPLE 


M X Tmeasured * 8 


1.0017233 
1.0012808 
1.0022497 
1.0012808 
1.0013437 
1.0019846 
1.0017433 
1.0017471 
1.0026321 
1.0019560 
1.0017424 
1.0026770 
1.0026207 
1.00131B0 
1.0026B46 
1.0024166 
1.0021400 
1.0025959 
1.0018959 
0.9916704 
1,0018024 
1,0019808 
1,0011015 
1.00133B0 
0.9986798 
1.0022001 
1.001B101 
1,0024242 
1.0010777 
1.0049429 
1.0026007 
1.0010719 
0.99993B0 
1.0010719 


-0.4258500 
-0.4281250 
20.4744628 
-0.4281250 
-0.2961913 
-0.399755B 
-0.2480956 
-0.3687988 
-0.3038086 
-0.3760253 
-0.9492187 
-1,0715818 
-0.9018085 
-0.8735839 
-0.8721191 
-0.8971679 
-0.6605957 
-0.8805175 
-0.6086425 
2.2940912 
-0.6160646 
-0,4994628 
-0.2665800 
-0.3956054 
0.3308593 
-0.4279296 
20.0564453 
-0.2967773 
-0.1212890 
-0.7456101 
-0.4801757 
-0.4208984 
-0.1928710 
-0.620B9B4 


THERMOCOUPLE 
THERMOCOUPLE 
THERMOCOUPLE 
THER MOCOUPLE 
THERMOCOUPLE 
THERMOCOUPLE 
THERMOCOUPLE 
THERMOCOUPLE 
THERMOCOUPLE 
THERMOCOUPLE 
THERMOCOUPLE 
THERMOCOUPLE 
THERMOCOUPLE 
THERMOCOUPLE 
THERMOCOUPLE 
THERMOCOUPLE 
THER MOCOUPLE 
THER MOCQUPLE 
THERMOCOUPLE 
THERMOCOUPLE 
THERMOCOUPLE 
THERMOCOUPLE 
THERMOCQUPLE 
THERMOCOUPLE 
THERMOCQUPLE 
THERMOCOUPLE 
THERMOCOUPLE 
THERMOCOUPLE 
THERMOCOUPLE 
THERMOCOUPLE 
THERMOCOUPLE 
THERMOCOUPLE 
THERMOCOUPLE 
THERMOCOUPLE 


1.0028639 
1.0017776 
1.0019627 
1.0019207 
1.0026245 
1.0024214 
1.0017195 
1,0013046 
1.0019274 
1.0026417 
0.9997830 
1.0019827 
1.0017757 
1.0030603 
1.0023384 
1.0024405 
1.0010166 
1.0015094 
1.0013409 
1.0013094 
1.0019588 
1.0011254 
1.0017185 
1.0252902 
1.00262B3 
1.0027475 

1.0016108 

1.0025406 
1.0019932 

1.0020733 

1.0018148 

1.00181«48 

1.001B339 

1.0009155 


-0.5545898 
-0.4930175 
-0.4770995 
-0.5316894 
-0.50332013 
-0.7976074 
-0.6256835 
-0.7516601 
-0.7330077 
-0.9266112 
-0.2326660 
-0.8005417 
-0.6931152 
-0.9011718 
-0.7065917 
-0.6196777 
-0.3536133 
>0.5520507 
-0.4960957 
-0.5520507 
-0.5768066 
-0.5678222 
-0.5253417 
-0.8309081 
-0.7043456 
-0.6196288 
-0,6912597 
-0.7730468 
-0.5812500 
-0.6718749 
-0.5980468 
-0.5980468 
-0.7208007 
-0.5206543 


Two plots of T actual vs Tmeasures for two thermocouples are 


presented in Figure A.1 and Figure A.2 to demonstrate their 


linearity. 
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APPENDIX B 
TABULATED DATA 


Heated Cylinder in the First Row 


По 211109 1.5m/s 3.0m/s 
Downstream Row No. Theta Theta 
1 1.0 1.0 
2 .468 „524 
3 23 46 .380 
4 2251] .311 
5 2191 1211 


Heated cylinder in the first row. 


Н = 4.6 1.5m/s 3.0m/s 
Downstream Row No. Theta Theta 
1 10 10 
2 .497 .484 
3 2365 .380 
4 276 277 
5 197 5 


Heated cylinder in the second row. 


Н = 2.3 1.5m/s 3.0m/s 
Downstream Row No. Theta Theta 
1 120 0 
2 501 .514 
3 .324 .354 
4 250 269 


S 


TABULATED DATA 


Row Ts o(°C) T (°C) Q(w) hlw/m2°c) Rey Nu U (m/s ) 
No. 


1 96.5 20785 5.03 56.61 394.2 14.09 8 
2 98.9 20.3 3. 05 56.71 394.2 14.11 888 
3 99.6 205 28 61 16 393.9 15.2 1 INS 
4 97.6 20.0 s. 60.34 394.5 15.02 888 
6 96.3 20.1 3.74 73507 394.4 18.39 5 
1 101.2 20.8 4.59 8 227 787.2 21.699 e 
2 89.4 20:7 4.17 93.26 787.4 23.19 8 
3 78.9 20.5 о 96.52 787.9 24.51 я: 
4 96.8 200612 4.74 95:26 787.6 23.7116 0 
6 86.2 20.3 4.57 107.69 788.3 26.900] M 
| 111:0 200 4.12 45.77 907.4 26.21 Е 
2 105.7 20-1 4.04 47.93 907.1 27.44 1.5 
3 112225 20.1 4.13 S105 907.1 29.22 S 
6 25 шиг! 3.78 52:61 907.1 30.10 Е 
1 99.3 2052 5.50 73.19 1813.7 41.90 M 
2 1071 20.1 5.99 69.86 1814.2 39.99 9 
3 112 20.2 6.81 73.20 1813.7 41080200 
6 а 201 oe ll 80.03 1814.2 45.9218 
1 25 20.0 3.27 42.74 1814.7 48.96 1.5 
2 QM 20:12 32027 42.49 1813.7 48.05 Е 
3 104.8 20.3 5.71 44.20 1813.2 50.6005 
4 100.7 20.6 3209 369 1811.6 42.22 888: 
6 Олио 20720 292 41.14 1811.6 47.06 1.5 
1 12122: 195 5.60 62.83 3634.7 72.040 D 
2 996 19.8 5402 65.70 3631.6 752528 IN 
3 956 12755 4.80 66.23 3630.6 75.70 008) 
4 20111572 20-1 3.92 56.91 3628.5 65.17 888 
6 85.8 20.1 9092 61.93 3628.5 70.91 И 
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